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Abstract 

We consider the Diophantine problem of Frobenius for semigroup S (d 3 ) where d 3 denotes the 
tuple [dx,d2,ds), gcd(di, cfe, cfo) = 1- Based on the Hadamard product of analytic functions [17] 
we have found the analytic representation for the diagonal elements atk (d 3 ) of the Johnson's 
matrix of minimal relations [12] in terms of di, cfo, g?3- Bearing in mind the results of the recent 
paper [10] this gives the analytic representation for the Frobenius number F (d 3 ), genus G(d 3 ) 
and the Hilbert series H(d 3 ; z) for the semigroups S (d 3 ). This representation does complement 
the Curtis' theorem [13] on the non-algebraic representation of the Frobenius number F (d 3 ). 
We also give a procedure to calculate the diagonal and off-diagonal elements of the Johnson's 
matrix. 
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1 Introduction 

Let S (di, . . . ,d m ) C N be the numerical semigroup generated by a minimal set 1 of integers 
{di, . . . ,d m } such that 

m < di < . . . < d m , gcd(di, . . . , d m ) = 1 , m > 2 . (1) 

For short we denote the tuple (di, . . . ,d m ) by d m . The least positive integer (d±) belonging to 
S (d m ) is called the multiplicity. The number (m) of the minimal generators of S (d\, . . . ,d m ) 
is called the embedding dimension. The conductor c (d m ) of S (d m ) is defined by c (d m ) := 
min {s € S (d m ) |s + NU{0}cS (d m )}. The genus G (d m ) of S (d m ) is defined as the cardinal- 
ity (#) of its complement A in N, i.e. A (d m ) = N \ S (d m ) and 

G(d m ) := #A(d m ) . (2) 

Introduce the generating function $ (d m ; z) for the set A(d m ) of unrepresentable integers 

$(d m ;z)= Z °- ( 3 ) 

s e A(d m ) 

The semigroup ring k[X\, . . . ,X m ] over a field k of characteristic associated with S (d m ) is a 
polynomial ring graded by degXi = di, 1 < i < m and generated by all monomials z di . The Hilbert 
series H(d m ; z) of a graded ring k [z dl , . . . , z dm ] is defined by [2] 

where Q(d m ; z) is a polynomial in z. Thus, we have 

H(d m ;z) + ^d m ;z) = T ^-, z<l. (5) 

The number F (d m ) := —1 + c(d m ) is referred to as Frobenius number. As follows from (3) 

F(d m ) =deg$(d m ;z) , G(d m )= ^ l s = $(d m ;l). (6) 

s e A(d m ) 

The determination of F(d m ), G (d m ) and i?(d m ;z) is called the m-dimensional (mD) Frobenius 
problem. The last restriction (m > 2) in (1) is essential since one can show that in the one 
dimensional case the problem is trivial 

H(d; z) = Y~^d ' F ^ = G ^ = 00 ' d - 2 • ( 7 ) 
The first non-trivial case (m=2) was studied already by J. Sylvester [3] 

H ^ Z ) = (i- z rf l)( i- z rf 2 ) ' F ( d ') = d ^ - d ±- d 2> G ( d2 ) = 2 (dl " 1)(d2 " X) • (8) 

The next non-trivial case (m=3) was extensively studied in the contents of commutative algebra 
[4]-[8] and algebraic geometry of monomial curves [9] where a partial progress was achieved (without 

lr The set {di, . . . ,d m } is called minimal if there are no nonnegative integers 6 4j for which the following linear 
dependence holds di = Yl^i bi,jdj, bij € {0, 1, . . .} for any i < m. It is classically known [1] that di > m. 
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calculating the Hilbert series). A new diagrammatic procedure of construction of the set A(d 3 ) was 
developed recently in [10]. It has paved the way to calculate Q (d 3 ;z) and, in accordance with (5) 
and (6), led to the complete solution of the 3D Frobenius problem. Based on Brauer's lemma [11] 
on the matrix representation of the set A(d 2 ) and Johnson's theorem [12] on the minimal relations, 
the author [10] was able to find the Frobenius number F (d 3 ) , genus G (d 3 ) and the Hilbert series 
H{d?; z) of a graded ring for the non-symmetric and symmetric semigroups S (d 3 ) : 

3 / 3 3 \ 

F(d 3 ) = I[(a,d) + J(d 3 )] G(d 3 ) = Ml + (a,d)-^4-I]a tt , (9) 

i=k V k=l k=l J 

Q n (d 3 ; z) = + Z 1 /2[(a,d)-j( d 3)] + z l/2[(a,d>+j(d3)] (1Q) 



k=l 



J 2 (d 3 ) = (a, d) 2 - 4 ^2 akkdjjdkdj + Uid 2 d 3 , (a, d) = ^ a kk d k , (11) 

i>j k=l 

Q s (d 3 ;z) = (\-z a ^) (l- z a ^) , if a 11 d 1 = a 22 d 2 , (12) 

where subscripts "n" in (10) and "s" in (12) stand for non-symmetric and symmetric semigroups 
[10], respectively. The matrix (((%•)) was introduced by Johnson [12]. Its elements uniquely define 
the minimal relations for given d\,d 2 , d 3 

audi = a\ 2 d 2 + ai 3 d 3 , 022^2 = 021^1 + 023^3 , 033^3 = a 3 i<ii + a 32 d 2 , where (13) 
an = min {tin I v u > 2, vudi = v 12 d 2 + ^13^3, ^12,^13 € N U {0}} , 
a 22 = min{i>22 I ^22 > 2, ^22^2 = v 2 \d\ +""23^3, v 2 i,v 23 G NU {0}} , (14) 
a 33 = min {v 33 I v 33 > 2, ^33^3 = v 31 di + v 32 d 2 , v 3 i,v 32 <E N U {0}} , and 
gcd(an,ai2,ai 3 ) = 1 , gcd(o 2 i, a 22 , a 23 ) = 1 , gcd(a 3 i, a 32 , a 33 ) = 1 . 

Notice that in formulas (9) - (12) only the diagonal elements of the Johnson's matrix (14) appear. 
Also notice that formulas (9) - (12) contain algebraic functions of ajj and dj, 1 < j < 3. However 
the ajj cannot be algebraic functions of dj because of the following theorem of Curtis [13] : 
There is no finite set of polynomials with integer coefficients 

{fj{xi, • • -,x m )} , j = 1, . . . ,n 

such that for each choice of d m ,m > 3 ; there is j such that fj(d±, . . . ,d m ) = F(d m ). 

On the other hand, it would be pretty interesting to build such representations for F (d 3 ) , G (d 3 ) 
and H(d 3 ; z) in the frameworks of analytic number theory avoiding any auxiliary invariants like the 
minimal relations (13) and (14). 

Our main result is the analytic representation for the diagonal elements a kk (d 3 ) of the Johnson's 
matrix (14) in terms of di, d 2 and d 3 : 

3 1 . ln* fc (d 3 ;z) 3 1 (1 - e id i dlt ) dt 

akk (d > ~ 4 hT^ ' ,Z >~2^J Q (l_ e ^)(l_ e ^)(l-^e-^) _1 ' 

(15) 

where (k,j,l) = (1,2,3), (2,3,1), (3,1,2). 

By (9) - (12) this gives the Frobenius number .F(d 3 ), genus G (d 3 ) and the Hilbert series 
H(d 3 ; z), and due to (5) it leads to the generating function <E> (d 3 ; z) for both symmetric and non- 
symmetric semigroups S (d 3 ) . 
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2 Matrix representation of the set A (d 2 ) and the map r 

First, recall the main statements about matrix representation of the set A (d 2 ). 

Lemma 1 ([11]) Let d\ and d 2 be relatively prime positive integers. Then every positive integer s 
not divisible by d\ or by d 2 is representable either in the form s = xd\ + yd 2 , x > 0, y > or in the 
form s = d\d 2 — pd\ — qd2, p, q £ N. 



Definition 1 Let 2 < d\ < d 2 with gcd(<ii, d 2 ) = 1. Define function a(p, q) as follows 

a(p,q) :=d 1 d 2 -pd 1 -qd 2 , p,q£N. 



(16) 



Lemma 2 ([10]) Let t be an integer and d 2 > d\. Then t £ A (d±,d 2 ) iff t is uniquely representable 
by 



1 <P < 



do- 



d L 



t = o~(p,q) , where 



1 < q < d\ — 1 , and d\ — 1 < 



d 2 



di 



(17) 
(18) 



where [c\ denotes an integer part of the number c. 



Lemma 3 ([10]) Let 2 < d\ < d 2 and gcd(d\,d 2 ) = 1. Every integer d% € A(d 2 ) gives rise to the 
minimal generating set {d\,d 2 ,d^} for the semigroup S (d 3 ) such that F(d 3 ) < F(d 2 ). 

The representation (17) of all integers a(p,q) <G A(d\,d 2 ) is called the matrix representation of the 
set A(di,d 2 ) and is denoted by M {A(d\,d 2 )} (see Figure 1), so 



a{M{A(d 1 ,d 2 )}} = A(d 1 ,d 2 ) , 



(19) 



i.e. a maps the entries of M {A(di, d 2 )} onto A(di,ci2) in a bijective manner (Lemma 2). 

a(p, q) is the integer which occurs in the row p and the column q of M {A(di, d 2 )}, e.g. <r(l, 1) = 
d\d 2 — d\ — d 2 . Note that the latter coincides with the Frobenius number F(d\,d 2 ). 



i 

2 
3 



F 2 


F 2 ~d 2 


F 2 -2d 2 






I', <l, cl. 


i;-(i r 2(i. 








ls-2(l,-2cl. 





















cl , —2 (I, —1 






q 




2^-d^ 





—if 











^<l, 2d, 









d l< 2. 



Figure 1: Typical matrix representation M {A(d±, d 2 )} of the set A(d±, d 2 ). The lowest cells in every 
column of M {A(di, d 2 )} (gray color) is occupied exclusively by the integers 1, . . . ,d± — 1 not in a 
necessarily consecutive order [10]. 

Throughout the paper we will make use of another entity, a map r, which was introduced in [10] 
in order to constitute the relationship between the set A (d 3 ) and the generating function $ (d 3 ; z) . 
Actually, the map r serves for such relationships in any dimension m > 1 as well. 
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Definition 2 The function r maps each power series (polynomials including) 

£c 8 [A] z s eN[z] 
sen 

with c s [A] G {0, 1} onto the set A of degrees {s € A C N I c s [A] / 0}. 
The coefficient c s [A] is a characteristic function for the set A C N 

which satisfies the following properties : 

Let two sets, A x and A 2 , be given such that A 1? A 2 C N. Then for set intersections and set 
unions, we have [14] 

c s [Ai n A 2 ] = c s [Ai] c s [Ai] , c s [Ai U A 2 ] = c s [Ai] + c s [Ai] - c s [A x n A 2 ] . (21) 

It also follows from Definition 2 and from representations (3) and (4) that 

r [$ (d m ; z)] = A (d m ) , r [# (d m ; *)] = S (d m ) , r [0] = , (22) 

T- 1 [A{d m )} = $(d m ;z) , r- 1 [S(d m )] =H(d m ;z) , r" 1 [0] = , (23) 

where denotes an empty set. Observe that since all coefficients of the polynomial <E> (d m ; z) are 
1 or 0, we can uniquely reconstruct &(d m ;z) from A (d m ) and vice versa. In this sense r is an 
isomorphic map. The map r is also linear in the following sense : 

Lemma 4 Let d m be given and a set A (d m ) be related to its generating function $ (d m ; z) by the 
isomorphic map r defined in (22). Let A\ and A 2; be subsets A (d m ). Then the following holds 

t- 1 [Ax U A 2 ] = r- 1 [Ax] + r- 1 [A 2 ] - r- 1 [Ax n A 2 ] . (24) 

Proof By (3), (20), (21) and Definition 2 we have 

t-^AiUAz] = Yl c s [A 1 UA 2 ]z s = (cs[Ax] + c s [Ax]-c s [AinA 2 ])z s = 

seM SEN 

r- 1 [Ai] + r- 1 [A 2 ] - r- 1 [Ai n A 2 ] , 

that proves the Lemma. □ 

The following statement is a consequence of Lemma 4 and generalizes the result obtained in [10] 
(formula (64)) 

Lemma 5 Let d m be given and a set A (d m ) be related to its generating function <3? (d m ; z) by the 
map t according to (22) and (23). Let n sets Aj, i = 1, . . . ,n be given such that 

A iC A(d m ), Aj n Aj = , for i^j = l,...,n. (25) 

Then the following holds 



.i=i 



Er-MA,] . 



(26) 



i=i 



Proof By induction on n applying Lemma 4 to the left hand side of (26) consecutively and making 
use of (25) we obtain 



n-l 



+ r- 1 [A n ]-r- 1 



1=1 

n-l 

\J {Ai n A n } 



n— 1 



+ T- 1 [AJ-T- 1 



' n— 1 



U Ai r fl A ™ 



. i=l 



1=1 



n— 1 

.i=i 



+ r- 1 [A n ] = ... = E^ 1 [A 4 ] 



i=i 



that proves the Lemma. □ 
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2.1 The intersection set A (d±, d 2 ) fl S (d 3 ) and its generating function 



In this Section we construct the intersection set A (di, d 2 ) H S (d 3 ). This set is essential to determine 
an explicit non-algebraic expression for a 33 (di, d 2 , d 3 ) (see next Section). 
First, define a ID numerical semigroup S (d 3 ) of integers a = mod (d 3 ) 



S (d 3 ) := {a \ a = jd 3 , j G N U {0}, d 3 > 2} , 



(27) 



which is generated by the Hilbert series H (d 3 ; z) given by (7). 

Consider the intersection set A (di, d 2 ) nS (d 3 ) which consists exclusively of the integers a = jd 3 , 
the index j runs with jumps from 1 to iV 3 such that N 3 d 3 G A (di,d 2 ) still holds. Notice that iV 3 
satisfies 

AT . did 2 -di-d 2 

N 3 < < di - 1 , 

that follows from iV 3 d 3 < F(d\,d 2 ) and d 2 {d\ — 1) — d\ < d 3 (di — 1). 

The integers jd 3 are distributed inside the matrix representation M {A(di,d 2 )} as shown in 
Figure 2. Let us prove an important Lemma. 



3 q 



N 3 cl . 



1 1 1 1 



2cl 



— r 



3d. 



d 3 



rcl3 



33 



a 33 < r, 11 < TSr 3 < tl , - 1 



Figure 2: Integers jd 3 G A (di,d 2 ) H S (d3) (black boxes) are distributed inside the matrix represen- 
tation M{A(di,d 2 )}. 



Lemma 6 Let d 3 &e given, d 3 = (di,d 2 ,d3) ; and assume that the minimal relations are defined by 
(13) and (U). Then 



and 



jd 3 G A (di, d 2 ) n S (d 3 ) , j = 1, . . . , a 33 - 1 , 



a 33 d 3 A(di,d 2 )nS(d 3 ) 



Notice that j runs in (28) without jumps. 

Proof The proof follows from two obvious relations. First, by (27) we have 

jd 3 G S(d 3 ) , j = l,...,a 33 ,... , 
and, next, by the 3rd line of (14) we have 

jd 3 G A (di, d 2 ) , j = 1, . . . , a 33 — 1 , and a 33 d 3 A (di, d 2 ) 



(28) 



(29) 



(30) 



(31) 
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that leads to (28) and (29). □ 

Keeping in mind Lemma 6 we write the generating function of the set A (d±, d 2 ) H S (g^). 
According to (23) we get 

033-1 ^3 

r- 1 [A (d 1; d 2 ) n S (d 3 )] = £ ^ + £ ^ > ( 32 ) 

j = l i>Q33 

where index j runs with jumps in the second summation of (32). Define an auxiliary function 
* 3 (d 3 ;z) by 

( d3 ' z ) = ~ T ~ l [A (dl> d2) n s • (33) 

It plays a key role in calculating the diagonal element 033. 

Theorem 1 Let d 3 be given, d 3 = {d\,d2,dz), and assume that the minimal relations are defined 
by (13) and (14)- Then the diagonal element 033 (d 3 ) is given by 

Md 3) = i_ U (34) 
v y d 3 z^o In z 

Proof A series expansion for ^3 (d 3 ;z) reads 

j<N 3 00 

$ 3 (d 3 ; z) = z a33ds + £ z jd3 + £ z id3 , (35) 

j>C33 j=Ar 3 +l 

where index j runs with jumps in the second summation of (35). The power series (35) converges 
for z < 1 

i<^3 JN 3 +i)d 3 
* 3 (d 3 ;z)=z^3 + ^ • ^-j- . 

i>«33 

The first term in (35) turns to be the leading term if z — ► 0, i.e. ^3 (d 3 ;z) ^— ^> ^^33^3 _ Hence, 
formula (34) follows. □ 

Formulas (33) and (34) have nothing specific concerned with arrangements of the di in the tuple 
(di,d 2 ,d 3 ). Therefore the similar formulas can be obtained for the diagonal elements an (d 3 ) and 
a 22 (d 3 ) 

, 1 ln* fc (d 3 ;z) , 

a tt d 3 =-lim \ K > , k = l,2, 36 

' dk In z 



z di z d2 



where 

* 1 (d 3 ;z)= T ^-r~ 1 [A(d 2 ,d 3 )nS(d 1 )] , ^ 2 (d 3 ;z) = -r- 1 [A(d 3 ,d 1 )nS(d 2 )] . 

3 Analytic representation of generating functions 

In the previous Section we have found the diagonal elements ajj (d 3 ) , j = 1,2,3, through the 
auxiliary functions tyj (d 3 ;z) ,j = 1,2,3. Thus, the non-algebraic expressions for ajj (d 3 ) will be 
given in the closed form if the analytic representation of corresponding generating functions 

T-MA^dOnStdfc)] , (k,j,l) = (1,2,3), (2,3,1), (3,1,2), 

will be found. The present Section is completely devoted to this question. In order to resolve it we 
make use of the Hadamard product of analytic functions. 
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3.1 Hadamard product of analytic functions 

Let u(z) = Y^=o a n,z n and v(z) = YlnLo bnZ n be analytic functions in the unit disk \z\ < 1. Define 
the <8>-product of u(z) and v(z) by 

1 t-2-K 00 

u(z)®v{z) = — u(ze it )v(e- it )dt = y" a n b n z n , \z\ < 1 . (37) 
27rJ ° ^0 

This product was introduced by J. Hadamard [15] to discuss the singularities of the analytic function 
f(z) defined by the series Y^=o a nb n z n in terms of those of the functions u and v ("multiplication 
of singularities" theorem [16]). 

Another representation, which is similar to (37), can be found for the analytic function f(z 2 ) = 
Y^=o a nb n z 2n through the integral convolution "o" of the functions u and v, 

1 f-2TT 00 

u(z) o v(z) = — u (ze u ) v (ze- u ) dt = Y* a n b n z 2n , \z\ < 1 . (38) 

This product was also used by J. Hadamard [15] and is strongly related to the <g>-product 

u (z 2 ) ® v (z 2 ) = u(z) o v{z) , (39) 

e.g. for the analytic function u(z) we have 

u(z) <g> (1 - z)~ l = u(z) , but u(z) o (1 - z)- 1 = u (z 2 ) . (40) 

Both (gi- and o-products have been variously referred to as the Hadamard product, quasi inner 
product or the termwise product (see survey [17]). The (8>- and o-products appear naturally in 
functional analysis, e.g. the Bieberbach conjecture for univalent analytic functions [18] and the 
Polya-Schoenberg conjecture for analytic convex mappings [19]. 

There is a variety of algebraic and other analytic properties associated with o- and 0-products 
[20]. The commutative and distributive laws can be easily verified 

u(z) o v(z) = v(z) o u(z) , u(z) o \v(z) + w(z)] = u(z) o v(z) + u(z) o w(z) , 

u(z) <8> v (z) = v(z) <S> u(z) , u(z) [v(z) + w(z)] = u(z) v(z) + u(z) <S> w(z) . (41) 

Notice that (38) manifests the non-associativity of the o-product 

[u(z) o v(z)\ o w(z) ^ u(z) o [v(z)] o w(z)\ , 

e.g. [z o z\ o z 2 = z 4 , but z o [z o z 2 ] = 0. However the 0-product is associative, 

[u(z) <g> v(z)] (8) w(z) = u(z) [v{z)\ (8) w(z)] , (42) 

e.g. [z <S> z] <S> z 2 =z®[z0z 2 ] = 0. The last property is very important to generalize the Hadamard 
product. We give an appropriate extension of 0-product referred to as the Hadamard multiple 
product [21]. 

Let Uk(z) = Y^=$ a-ri^z n , k = 1, . . . , N, be analytic functions in the unit disk \z\ < 1. Define the 
Hadamard multiple product of Uk(z) by 



N 



u k {z) = m(z) u 2 {z) <g> . . . <g> u N {z) = ^ a n ] -a {2) ■ a { n N) • z n . (43) 



k=l n=Q 



8 



The (^-product is well defined due to the associativity of the 0-product (42). Its integral repre- 
sentation has the form 

N 1 f 27T f 27T / N \ N 

= j^yjZT J o ■■■J o u^ze^^llukie-^dtk. (44) 
We give some examples of the <g>-product operating with rational functions [16] 

N „ N 

, n k € N , L = lcm (m, . . . , un) (45) 



and also of the o-product operating with transcendental functions [20], [21] 

z z z z z z z z 

cosh - o cosh - + sinh - o sinh - = Iq{z) , cosh - o cosh - — sinh - o sinh - = Jo(z) , (46) 

£ £1 £i £1 £ Zi £j £t 

e z/2 Q e z/2 = Iq{z) ^ e - g Q (1 + z) n = Ln(z 2 } j _L_ Q _±_ = ^ ( Q; . 1; ^ (4?) 

where a,f3 > 0. In (46) and (47) the functions Io(z), Jq(z), L n {z) and 2F1 {pt; (3; 1; z) denote 
the modified and non-modified Bessel functions of the 1st kind, the Laguerre polynomial and the 
hypergeometric function, respectively. 

Because of its connection with the integral convolution (38), the Hadamard product leads to 
elegant evaluations of complicated trigonometric integrals and provide analytic derivations of com- 
binatorial identities [20], [21]. Finally, the Hadamard product was used in [22] to construct solutions 
for a variety of Cauchy-type problems. 

It turns out that the <g)-product is a sensitive tool to deal with an intersection of discrete 
numerical sets. This will be the subject for discussion in the next Section 3.2. In order to be 
technically equiped we derive here two important formulas which will be useful further. 

Lemma 7 Let u(z) = YH^o^^ be an analytic function in the unit disk \z\ < |c|. Then 

u{z) ® — *— = - u (-) , c + . (48) 
c — z c \cJ 

Proof Calculating the <g>-product we get 

. . 1 1 . , 1 1 -r-^ u sr^( z \ k 1 z 1 f z \ 

u(z) ® = - u(z) - = - \a k z <g> > I - I =-> a k -r = -u[- \ , 

c — z c 1 — z/c c \ct c c K c \c/ 

' k=0 k=0 k=0 

that proves formula (48). □ 

Before going to the next formula (Lemma 8) we prove an elementary identity 

n ST^ w n (2m\ 

—1 = ^7^' ^ = ex HvJ • (49) 



Rewrite the right hand side of (49) in the form 

< _ 1 f- f ^k-iurr ^- fc _n 1 ^ 1 -2n 2 ^- 2 + ... + n(-i)«- 1 n r , 

=1 - W n z 1 fc=1 

where Yl k denote the basic symmetric polynomials 



f— ' z — wi z 11 — 1 f— ' z n — 1 



U k = W n W n ■ ■ ■ W n , ™- n i = E W ™ ^ = E W n W n > ' H - = II W n' ( 51 ) 

h>h>— >ifc=l 1=1 h>h=l 1=1 
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Recall Vieta's formula for the sum Sk of the products of distinct roots of the polynomial equation 
of degree n 

a n z n + a n -iz n ~ 1 + . . . + a\z + a = , 

which reads [23] 

n 

S k = (-1)*^ > wh ere S k = £ z,^ . . . z h . (52) 

™ «l>«2>--->ifc = l 

Since the u^-roots are associated with the polynomial equation z n — 1 = then the Vieta's formula 
(52) gives for Hfc -polynomials 

f , if 1 < k < n , 
(-1)" +1 , if & = n. (53j 

Substituting (53) into (50) we arrive at (49). 

Lemma 8 Lei u(z) = YlT=o a k zk be an analytic function in the unit disk \z\ < 1. XTien 

n-l 



1 1 

u ^®r3^r = - E n H0 • ( 54 ) 

fc=0 

Proof Making use of the identity (49) we obtain by Lemma 7 

fc=l V n 7 Z=l n V n/ Z=l fc=0 



that proves formula (54). □ 

Henceforth we use the ^-product and refer to as the Hadamard product. 

3.2 Representation of r" 1 [A (d 1} d 2 ) n S (d 3 )] 

In this Section we derive the analytic representation of the generating function for the intersection 
of two sets, A (di, 0I2) and S (cfe). Start with the following proposition which utilizes the Hadamard 
product for the intersecting sets Ai and A2. 

Lemma 9 Let Ai, A2 be subsets ofN, and let their corresponding generating functions r _1 [Ai] 
and r _1 [A2] be 

r- 1 [A l ]=^2c s [A i ]z s , i = l,2, (55) 

where c s [Aj] stands for corresponding characteristic function of the set Aj and satisfies (20). Then 
the set Ai PI A2 is generated by 

t- 1 [Ax n A 2 ] = r- 1 [Ax] ® r- 1 [A 2 ] . (56) 

Proof By Definition 2 of generating function r _1 [A] and the 1st relation in (21) we have 

r- 1 [Ax n A 2 ] = c * t A i n A 2] ^ = E c * t Al l c * [ Al l • ( 57 ) 
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On the other hand, by definition (37) of the Hadamard product we have 

r- 1 [Ax] ® r- 1 [A 2 ] = ^c s [Ax] z s ® J] c s [A 2 ] z fl = [Ax] c s [Ax] z s . (58) 

seN seN seN 

A comparison of (57) and (58) proves the Lemma. □ 
Lemma 9 has a simple generalization 

N 

T- 1 [n? =1 A j ]=<g)T- 1 [A j ] . (59) 

3=1 

Returning to t~ 1 [A (dx, d 2 ) H S (d 3 )] we can verify that both generating functions r _1 [A (dx, d 2 )] 
and r _1 [S (da)] are representable in the form (55) 

r- 1 [A(dx,d 2 )] = $({di,d 2 };z)= £ *«, T -i [ S (cfe)] = = £ z s . 

s e AfiiiA) ' ' s e S(d 3 ) 

Therefore we obtain 

r- 1 [A (d^ck) n S (d 3 )] = r- 1 [A (dx, da)] ® r" 1 [S (d 3 )] = $ ({dx, d 2 }; z) ® — . (60) 

Formula (60) can be slightly simplified by utilizing the relationship (5) between the generating 
function $ and the Hilbert series H and of distributive law (41) for the Hadamard product 

r- 1 [A (dx, d 2 ) n S (d 3 )] = ^ ® 7^7 " ^12 W ® 73^ > (61) 

where H\2(z) denotes for short the Hilbert series H ({dx, d 2 }; z). On the last step we make use of 
the 1st equality in (40) and get the generating function for the intersection of two sets, A (dx,d 2 ) 
and S (d 3 ), 

r- 1 [A (dx, d 2 ) n S (d 3 )] = - ff 12 (s) ® — L- , (62) 

and also the function vl/3 (d 3 ; z) introduced in (33) 

* 3 ( d3 ; z ) = r-^h - T ~ l t A d ^ n s ^ = H ^ ® - 1 • ( 63 ) 

Applying now Lemma 8 to formula (63) we get 

*3 (d 3 ; z) = 1 £ i?i2 (z<) - 1 , ^3 = exp (^) . (64) 
Simplifying the expression (8) for H12 (z) 

p=0 p=0 g=0 

and substituting it into (64) we obtain 

^ d'2 — 1 OO (^3 — 1 

* 3 ( d 3. z ) = J_ ^ ^ ^1+^2 £ w Kpd 1+qd2 ) _ 1 (g6) 
3 p =o 9 =0 fc=0 
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The inner sum in (66) does vanish for p, q such that pd\ + qd 2 is not divisible by dj,. Indeed, 

V- 1 k(pd 1+ qd 2 ) exp [2m (pdi + qd 2 )] - 1 „ 
fc=o exp (27ri £ ^ J - 1 

Thus, it remains 

d2 — 1 OO — 1 CO oo 

^ 3 (d 3 ;z) = E z pdl+qd2 ~ 1 = Yl z pdl+qd2 =J2 zbjd3 , ( 67 ) 

p=0 9=0 p=l 9 = 1 j=l 

d ^\P d l+1 d 2 <%|j>d 1 +gd 2 

where 6j G N is defined as the integer which satisfies the Diophantine equations in bj 

pd\ + qd2 = bjds , p = 0, . . . , cfe - 1 , <Z = 0, . . . (68) 

at least with one solution. Recalling the definition (13) and (14) of the Johnson's minimal relations 
we conclude that 

bi = a 33 , (69) 

and the first term in series expansion (67) reads z a:i:id:i . Hence, formula (34) follows. Notice that 
(67) can be obtained by straightforward calculation of the Hadamard product (63) according to 
Lemma 9 and representation (65) 

(d 2 — 1 oo \ oo /cfc — 1 oo \ oo oo 

Y E zPdl+qd2 ® E z '' &i - 1= EE zPdl+qd2 ® E zTdA = E zbjd3 ' 
p=0 q=0 J r=0 \P=1 9=1 / r=l j=l 

where bj is defined in (68). 

We finish this Section by giving an integral representation for (63) according to (37) 

3 1 f 2n (1 - e idld2t ) dt 

^ 3 ^ ' Z > = 27 J (1 - e^') (1 - e*^) (l _ z d ie -id z t) ~ 1 • ( 70 ) 

The other two functions, (d 3 ; z) and ^2 (d 3 ; z) , can be written by the cyclic permutation of the 
indices (1,2,3) in (70) 

3 1 f 2 * (1 - e id2d3t ) dt 

* 1 ^ ' Z > ~ 27 J (1 - e id2t ) (1 - e^3*) (l _ z d, e -id lt ) ~ 1 ' ( 71 ) 

3 1 f 2n (l - e id3dl *) dt 

^ 2 ^ d ' ^ = 27 / (1 - e id ^) (1 - e id i') (1 - ^e-ifct) " 1 ' (72) 

4 The explicit calculation of the entries of the Johnson's matrix 

The complexity of expressions (34) and (36) in conjunction with (70) - (72) makes further evaluation 
of the diagonal elements akk (d 3 ) of the Johnson's matrix excessively difficult. Therefore, we develop 
in this Section another approach representing akk (d 3 ) as zeroes of some functions. 
Combining (67) and (70) we have 



oc 



j_ r (i ~ eidid2t) dt i=v.M3 m) 

2vr J (1 - e id ^) (1 - e^*) (1 - ^ 3e -id 3 t) ^ 
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where b±, b%, . . . present all positive integers such that in accordance with (68) every integer bjd^ is 
representable by d\ and d.2- 

Fix the index j = k and differentiate with respect to z ds both sides of the equality (73) bk times, 
and take their values at z = 



b_l ^ (1 

2vr 



id\d2t 



) 



-ib^dzt 







-dt = h\ 



Define a new function S3 (d 3 ; 6) by 



1 



]3. 



2- 



) 



-ibd^t 



(\ — gicM^ M _ e id 2 t 



-dt . 



(74) 



(75) 



(76) 



Thus b k is its zero, S3 (d ; 6^) = 0, and 033 (d ) is its minimal integer zero according to (69) 

a 33 (d 3 ) = min {b k \ E 3 (d 3 ; b k ) = 0, b k G N} . 

Notice that 2 < 033 (d 3 ) < d\ — 1 according to [10], that bound the range of b where the first zero 61 
does appear. In Figure 3 we present the typical plot of the function S3 (d 3 ; b) for the triple d\ = 23, 
c?2 = 29 and d% = 44 which was considered numerically in [24] 





1 
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. 6 





. 4 





. 2 












10 
to 



• •••••• 9^ 



15 
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Figure 3: Typical plot of the function S3 ({23, 29, 44}; b) and the distribution of its zeroes b k . The 
diagonal element 033 of the Johnson's matrix reads 033 (23, 29, 44) = 5. 

The other two diagonal elements, an (d 3 ) and 022 (d 3 ) , can be obtained by the cyclic permuta- 
tion of the indices (1,2,3) in (75) and (76) 

a n (d 3 ) = min{o fe I Si (d 3 ;o fc ) = 0, b k G N} , a 22 (d 3 ) = min{6 fc I S 2 (d 3 ;6 fc ) = 0, b k G N} (77) 
where 

1 f 27T (l — e id2dsts ) e~ ibdlt 3 1 f 2n (l — e id3dlt ^ e~ ibd2t 

2^ J (1 - e id2t ) (1 - e id ^) dt ' "2 (d , 6) — 1 - — j^—^^—^-^dt . 

and 2 < a u (d 3 ) < d 2 - 1 , 2 < a 22 (d 3 ) < di - 1. 

Performing the calculation for the triple di = 23, ^2 = 29 and d3 = 44 we can find in accordance 
with (77) 

an (23, 29, 44) = a 22 (23, 29, 44) = 7 . 

Together with 033 (23, 29, 44) = 5 these are in full agreement with the Johnson's matrix of minimal 
relations which was found in [10], Example 2. 
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4.1 The off-digonal elements of the Johnson's matrix 

The uniqueness of the Johnson's matrix ((a^)) of minimal relations makes us possible to determine 
also its off-digonal elements for non-symmetric semigroup S (d 3 ). As was shown by Johnson [12], 
six off-digonal elements of the matrix (((% )) are related by six identities 

021 + 031 = an , ai2 + a 32 = a 2 2 , 013 + 023 = 033 , (78) 

«23 a 32 = 022033 — ^1 j a 13 a 31 = 011^33 ~~ d 2 , «12«21 = Ona 22 — ^3 • 

These identities give rise to six quadratic equations 

^3^23 ~ (( a > d ) ~ 2andi) a 23 + (022033 - di) d 2 = d 3 a\ 3 - ((a, d) - 2a 22 d 2 ) a 13 + (a u a 33 - d 2 ) di = , 
^32 ~ (( a > d ) ~ 2a n rf i) 032 + (033022 - di) d 3 = d 2 a\ 2 - ((a, d) - 2a 33 d 3 ) a 12 + (ana 2 2 - d 3 ) d 1 = , 
di a ii ~ (( a > d ) ~ 2a 22d 2 ) 031 + (a 33 a u - d 2 ) d 3 = dxa 2 21 - ((a, d) - 2a 33 d 3 ) a 21 + (a 22 an - d 3 ) d 2 = , 

where (a, d) is already defined in (11). Notice that all these equations have common discriminant 

(a, d) 2 - 4 (ana 22 did 2 + a 22 a 33 d 2 d 3 + 033011^3^1) + 4did 2 d 3 , 

which can be recognized as J 2 (d 3 ) defined in (11). As was shown in [10] its square root, J (d 3 ), is 
a positive integer, 

J (d 3 ) = |ai 2 a 23 a 3 i - ai 3 a 32 a 2 i| > 1 . (79) 
Therefore the solutions of all six quadratic equations are always rational numbers. 

4 = ^ ((a, d) ± J (d 3 ) - 2ondi) , 4 = 2^ « a > d ) ± J ( d3 ) " 2 «n d i) > 



1 . - -x 1 



°31 = 



((a, d) ± J (d 3 ) - 2a 22 d 2 ) , af 3 = — ((a, d) ± J (d 3 ) - 2a 22 d 2 ) 



2di vv ' ' v ' ' 1 16 2d 3 

af 2 = ((a, d) ± J (d 3 ) - 2a 33 d 3 ) , 4 = ^ « a , d ) ± ^ ( d ') " 2o 3 3d 3 ) • (80) 

Show that one of two rational roots associated with every quadratic equation is always positive 
integer while another is necessarily not integer. Making use of the identity [10], formula (135), 

- [(a,d) ± J (d 3 )] = ana 22 a 33 + - (ai 2 a 23 a 3 i + ai 3 a 32 a 2 i ± |ai 2 a 23 a 3 i - ai 3 a 32 a 2 i|) , 

and the identities (78), we have 

0^3 = °23 , a 3 2 = a 32 , a^i = a 3 i , af 3 = a 13 , af 2 = a 12 , = a 21 , (81) 

d 2 d 3 d 3 di di d 2 

a 23 — a 32 — , a 32 — a 23 — , a 31 — 013-7- > a 13 — a 3 i— , a 12 — a 2 i— , a 21 — oi 2 — . 
03 ~ (l> <l\ a 3 a 2 " 0. 



There is only one way to satisfy the uniqueness of the Johnson's matrix ((a^)) of minimal relations 



comprised exclusively of integer entries a™ if we require that a {j = af, = a™. However, as one can 



see from (80), this leads to J (d 3 ) = that contradicts (79). 
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